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ABSTRACT 

Satellite measurements of sea surface temperature must be corrected for 
atmospheric moisture, cloud contamination, reflected solar radiation and other 
sources of error. Procedures for reducing errors are discussed. It appears 
that routine accuracies of 1°C are possible, given low noise spectral measure- 
ments in the infrared. 
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I. Introduction 


Although present satellite measurements are capable of providing accurate 
estimates of sea surface temperature in local areas (Smith, et al. 1970; Rao and 
Smith, 1972), the achievement of routine 1°C accuracies on a world- wide basis - 
a specific global atmospheric research program (GARP) requirement — has not 
yet been demonstrated. This paper investigates some of the problems which 
tend to interfere with accurate determination of sea surface temperature (SST). 

It is likely that with proper satellite instrument design and data handling pro- 
cedures these error sources can be reduced to the desired levels. 

Section II deals with techniques for measuring the atmospheric correction 
to be applied to "window" radiance measurements, and the necessary instrument 
specifications to hold the effect of random system errors to a low level. 

Section III deals with the procedure suggested by Smith and Rao (1971) for 
using 3.7 and 11 micrometer spectral information to eliminate the effects of 
clouds. Recommendations for improvement are made, although the method still 
will not always produce unique answers. Section IV analyses the stability of sta- 
tistical regression coefficients as an indicator of accuracy for prediction. 

It must be noted that several types of error are not dealt with here. The 
first are system errors arising in calibration of the satellite instrument, 
telemetry, analog to digital conversion, satellite pointing errors, etc. For an 
analysis of a current system see the review by Leese et al. (1971). The second 
type of error arises from identifying radiance measurements with the nominal 
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"bucket" temperature used by oceanographers. The difference between the 
"skin" radiance temperature and the internal temperature may go over 1°C 
(Paulson and Parker, 1972), but it seems unlikely that this will contribute 
significant error on a global scale. 

Finally, reflected solar energy can cause a noticeable error (Whitehead, 
1972; Greaves, 1972) in ll^m radiances under very calm surface conditions. 

For this reason we recommend the rejection of data which come from within 
5° of the point of specular reflection of the sun. 

II. A. The Use of Two Infrared Measurements for Eliminating the Effect of 
Atmospheric Water Vapor 

Although measurements in the 11 micrometer atmosphere "window" yield 
reasonable estimates of SST, a correction for attenuation from atmospheric 
moisture must be applied if high accuracies are desired. The calculations 
described here show that the correction may be as large as 8°C for a very warm, 
moist atmosphere. Of course atmospheric water vapor is highly correlated with 
the temperature of the sea surface, so that the error in using 11/^m radiances 
arises from the application of a standard correction as a function of SST. Thus 
(Temperature of sea surface) = (Temperature from 11/^m radiance) 

+ AT(1 1/xm temperature) 

where AT(T) is a prescribed function. In practice the quantity AT is not 
determined solely by the surface temperature, but also by the prevailing 
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climate of the region, and by variations in the weather. Whether the error com- 
mitted by taking AT in (1) to be a fixed function amounts to 1°C on a global basis 
is not known. Certainly the neglect of time and space variations will lead to 
systematic errors (e.g. dry high pressure areas will yield high estimates of 
SST), which will interfere with long range weather forecasting. 

As described by Anding and Kauth, (1970) satellite measurements in two 
spectral windows, one "clean" and one "dirty," can provide additional informa- 
tion to make measurements accurate to 1°C on a case by case basis. The choice 
of these spectral c hann els is still open, as several factors must be taken into 
account. 

In order to establish quantitative procedures for deriving sea temperature, 
we start with the equation of radiative transfer, in the usual notation, 

I, = B V (T S ) t/p s ) - 

where l v is the observed radiance at frequency v . We may rewrite 2 as 

K = B ,(T S ) T „(P.) + [1 - n/Ps)) < B v> 

where y , the weighted atmospheric emission, is defined by 

J B tT( p )] — &?/w - T v ( p s )] 


I 


p ” 3 t (p) 

B[T(p) 3 — I dp 

Bp 


( 2 ) 


3 



In the linear approximation of small absorption appropriate to atmospheric 
windows we have 1 - r v (p s ) = = k v u, where u is the path length and e is 

much less than one. Also 


dr' 

V 


dp v dp 


where q is the mixing ratio. In this approximation 



r Ps 

B[t(p)] q(p) dp 

0 


and we note that the atmospheric emission is independent of the absorption 
properties of the water vapor. Thus 

K = B v( T s)- e vKB v >-B y (T s )] 


In the spectral intervals we consider the moisture is at low levels, and we 
may expand the Planck function B about the surface temperature. For two 
spectral channels measuring radiance temperatures and T 2 we have 

= ^sst “ € i^a 

^2 = ^SST " e 2^A 

where T is the effective atmospheric temperature. We solve for T and, 

A oT 

following MacMillan (1971), define y = € x /(e - e 1 ) 


^sst _ Tj + “/(T j - t 2 ) 


( 4 ) 
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This theoretical relation is complicated in practice by several factors: 

(1) spectral regions with stronger molecular absorptions tend to have 
strong line absorption, which does not obey the assumptions we have made. 

(2) there is considerable evidence (Bignell, 1970) that the absorption co- 
efficient k v for water vapor depends on q, so that the linear approximation is 
satisfied for only very small water vapor amounts. 

It appears that both of these effects can be accounted for simply by 
generalizing (4) to 

Tsst = T ! + ^ T i - T 2 ) + c (5) 

where C is a constant to be determined from comparison of satellite data with 
measured SST’s. 

Table 1 presents the results of calculations of the atmospheric correction 
for clear sky atmospheres taken from the study of Wark et al. (1962). The 
calculations were carried out using the high spectral resolution radiative trans- 
fer program developed by Kunde and Maguire (1973). The high resolution results 
were averaged across the spectral band to produce the figures indicated. 

The spectral interval 870 - 950 cm -1 (10.5 - 11.5 micrometers) has very 
little atmospheric absorption and has been used frequently on satellite instru- 
ments as an indicator of surface or cloud temperatures. The interval 775 - 831 
cm’ 1 (12.0 - 12.9^m) has been recommended by Prabhakara et al. (1972) for 
estimating the atmospheric correction for SST. This recommendation is 
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especially valuable as it is based on the examination of satellite data (Nimbus 
IV IRIS). The interval 1150-1250 cm' 1 (8.0-8. 7^m) was suggested by Kunde as 
a candidate because it is affected mainly by water vapor, but has stronger ab- 
sorption than the 11 micrometer region of the spectrum. 

By using a least squares fitting procedure we obtain the following equations 
for the surface temperature 


T = 

SST 

3.67 T„ - 

■ 2.84 T 12 s + 45.6 

rms error = 0. 7°C 

(6a) 

^SST = 

S-35T,., 

- 3.90T n - 105.6 

= 1.6°C 

(6b) 

^SST = 

3 27T 8.4 

- 2.23T 12 5 - 8.22 

= 0. 3 °C 

(6b) 


It is clear that the third case is the best choice for the spectral channels 
according to the model. This is not surprising as the contrast in the water 
vapor absorption is strongest for this choice of spectral intervals. 

Limitations on computer time prevented the running of more model calcu- 
lations. Due to the uncertainty in knowledge of the absorption coefficients it 
would be desirable to test the predictions above with an experimental program. 

With confidence in the qualitative behavior of the theoretical model (Kunde 
et al., 1973) we find that spectral information can be used to correct for attenu- 
ation from atmospheric moisture. 
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B. The Effect of System Errors on Two Channel Measurements. 

In addition to the errors from imperfect information as to the atmospheric 
state we must consider errors in the observing instruments, data handling, data 
processing, etc. From equation 4 we may estimate the error 8 (SST) arising 
from errors ST and ST 2 in the respective spectral measurements. 

8 (SST) = [8Tf(l + y 2 ) + 8T 2 y 2 - 2y(l + y) ST^] 1/2 

The quantity ST X 8 t 2 may be assumed to be zero, except for bias errors; ob- 
viously we wish to have both biases of the same sign (both channels high, or both 
low). Ignoring 8T ST 2 , and assuming ST 2 ~8T 2 we have 

S (SST) = [8T*(1 + 2y + 2y 2 )] 1/2 ( 7 ) 

Retracing the logic we see that we wish to have T x - T 2 as large as possible, in 
order to make y as small as possible. However, when this quantity is large 
there must be strong absorption in one of the channels, so that the linear approxi- 
mation fails. The tradeoff is between a strongly absorbing spectral region, in 
order to minimize the effect of system noise, and a weakly absorbing region 
such that the linear absorption approximation is justified. 

From 6c and 7 we see that 'y '-o 2.25 according to the model calculations. 

Thus error in the measuring factor is amplified by a factor /I + 2y + 2 y 2 ^3.8 
so that low noise measurements are essential. 
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Recently Anding and Kauth (1972) have carried out calculations similar to 
those described here. They recommend spectral channels at 8.9 A<.m and 11.9 pan. 
Based on their published data we may obtain an estimate of y ~3.6, which yields 
an error amplification factor of about 5.5. 

The necessity for low noise measurements is clearly indicated by these 
results. We are optimistic that these requirements can be met with suitable 
instrument design. 

in. Elimination of Cloud Effects by Combining Spatial and Spectral Information 

Clouds represent a persistent problem in satellite meteorology. Total 
coverage by high clouds prevents the use of infrared data to determine SST in the 
affected region. A less obvious problem is the fact that partial cloudiness, or 
complete coverage by very low clouds may cause slightly lowered (and hence 
false) surface temperature results. The development of procedures and 
thresholds for rejecting these cloud contaminated results is a subject of 
continuing research. 

Smith and Rao (1971) have described a method which uses pairs of spectral 
measurements (the 3.7 and 11 micrometer atmospheric windows) from spatially 
nearby points. The method rests on the assumption that clouds have constant 
spectral properties in a set of measurements in a given locale. The assumption 
is quite reasonable, and one can probably account for variations in cloud prop- 
erties (height, thickness, etc.) by increasing the estimated magnitude of random 
error in the measuring system. 
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Granted this assumption, we may use the nonlinear temperature dependence 
of the Planck function to solve for the surface temperature. This results from 
the fact that clear and cloudy areas contribute different fractions to the total : 
energy in different spectral intervals. 

Smith and Rao suggest using numerous pairs of measurements in a given 
region to obtain a series of estimates of SST. The average result is then their 
best estimate of the true sea temperature. Because of non linear error propo- j 
gation it is preferable to perform averaging on the raw data, and use the smoothed 
data in extracting the result. 

Figure la shows the behavior of radiance temperatures in two spectral 
bands as a function of the percentage cloud cover, assuming that clouds radiate 
as black bodies. The variation in T 3 7 - T n results from the fact that radia- . 
tion varies as a higher power of temperature in the 3.7ptm window. 

Calculations (Curran 1972) show that clouds have a lower radiation tem- 
perature at 3.7 Mm than at 11. This effect may be modeled by assigning the 

cloud an emissivity less than unity at 3.7m, with the result shown in figure lb. , 

! 

We s hall assign the cloud an emissivity of .75 at 3.7 m m. In figure 2 radiance : 
temperatures are plotted for two surface temperatures 280°K and 290°K for a 
number of cloud temperatures, with T/< T surface for all cases. This figure 
suggests an algorithm for obtaining the surface temperature from satellite data. 

It will not be difficult to include information from a 9 or 12 Mm channel in order 
to correct for molecular attenuation. In this discussion we ignore molecular 
absorption. 
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A. System noise tends to thicken the lines on figure 2. As described in 
section IIB the error in T 3 7 - T n is magnified because it is a dif- 
ferenced quantity. Furthermore we expect greater variability in cloud 
properties as cloud amount increases. For these reasons we identify 
warmest temperature, T , from a set of measurements, and discard 

CO 

all data which are much colder, say T &j - T > 10°. If too few data pass 
this test no result can be obtained for surface temperature. 

B. In order to average out random errors we fit a parabola through the 
rem a ining data. The root(s) where T 3 7 = T n £ represent possible 
values for the surface temperature. There are now two possible diffi- 
culties. 

1. The slope d(T 3 7 - T u )/dT 3 7 of the parabola is near zero in the 
vicinity of the point T = T u . In this case errors in the measured 
values result in greatly magnified errors in the estimate of SST. 

For cases where this slope is less than 0.3 we reject any extrapolation 
and fall back to the one channel histogram method of Smith (1970). 

2. In figure 3 the case of 270° clouds may be identified with SST = 280° 
for small cloud amount and with SST = 290° for large cloud amount. 
This ambiguity may occur when the slope is positive and T g 7 - T 
is negative. In this case again we revert to the one channel histo- 
gram for an answer. 
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For data processing on an operational basis this case may generally 
be resolved by comparing the predicted value with those obtained 
earlier in time. It is highly desirable to inject new results into a 
r unning time average in order to reduce the effects of random 
errors and to produce values of SST in regions where cloudiness 
causes frequent data gaps. 

C. Because the extrapolation by a parabolic fit tends to magnify errors in 
the data we require that a reasonable fraction, say 10% of the measure- 
ments fall within 5° of the predicted temperature. It is apparent that 
this chain of logic does not always yield an answer. The method works 
better for high, cold clouds than for low, warm clouds, and it tends to 
fail as the total cloud cover increases. This result is very plausible 
on physical grounds. 

By using the improved procedure described here results for SST can be 
obtained even in the presence of considerable cloud cover. 

IV. Stability of Statistical Correlation Coefficients as an Indicator of 

Prediction Capability 

Users of satellite data frequently have the task of deriving results from an 
incomplete set of measurements. For example a knowledge of low level water 
vapor is required for an accurate measurement of SST. This quantity is not 
determined by the measurements of current radiometers. 
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In this case one may tie together by empirical relationships the satellite 
measurements and independent ground truth measurements. These relationships 
may then be used to extend satellite results to regions lacking observations, or 
to produce results on a time scale for which ground truth data are not available. 

It is desirable to know the probable error resulting from the extrapolation 
of such statistical relations into regions in time or space where independent 
information is not available. 

For present purposes a simplified treatment is adequate. Let the desired 
quantity T be related to variables A and B by 

T = aA + /SB ( 8 ) 

and assume that our instrumentation measures only A. We may take advantage 
of correlation between A and B by writing 


T = 7A + C 


(9) 


where 7 and C are obtained by fitting the measured values of A to independent 
measurements of T. We use < > to represent an average over measurements 
used in the statistical analysis. We obtain 


7 = 


a + 


P 


<ab) - <aXp> 
_<A 2 >-<A^ _ 


C = P 


<B> - A 


<AB) - <A><B> 
<A 2 > - <A> 2 
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which illustrates the dependence of 7 and C on the correlation between A and 
B. From (9) and its average we have 


T - [<A><T>-<AT>] (A _ <A>) + <T> 

<^> - < A > J 

so by identity 

y - <A> <X> - < AT > 

A 2 - <A> 2 
C = <T> - y <A> 

We may now estimate the error in T arising from a variation of the fitting 
coefficients 

(ST ) 2 = StJ = (A - <A» 2 (Sy ) 2 

In a recent paper Shenk and Salomonson (1972) obtained equations relating 
SST to 11 ft m radiances and other spectral measurements. 

SST = .836 T n + • • • 31 - 37°N 

SST = 1.17T U + ■ • • 37 - 43°N 

Comparing with equation (9) we have Sy ~ (1.17 - .84) = .33 and since 
T n - SST, we may estimate the variability of T n from the variation of the 
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derived values for SST. From their figure 10 we find the variation in each 
latitude band is about 4°C. Thus /(ST) 2 -(4°C) x .33 - 1.3°C. 

This is a conservative (large) estimate for the possible error, as some of 
the variation of y undoubtedly arises from climatic differences in the two lati- 
tude bands. If we consider the variation of y about the average value, we find 
/ 8T 2 ~ 0.7°C. 

When other error sources are included such as instrument noise, errors in 
the ground truth data, etc., it is clear that a direct measurement of SST is 
preferable. 

VI. Conclusion 

It appears feasible to obtain sea surface temperature from satellite meas- 
urements, with an accuracy of 1°C. The effect of high cirrus (Braun, 1971) can 
be minimized by using 6.7/xm water vapor measurements. Additional information 
from the microwave region of the spectrum will prove useful in overcoming the 
effects of clouds. Other cloud effects and system errors can largely be eliminated 
by the procedures described here. A suitable low noise multichannel instrument 
is planned for future satellite missions. 
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Table 1 


Sounding (1200 GMT) 

Surface 

Temperatures 


Radiance 

Temperatures 


12.0 - 12.9 ^m 

10.5 - 11.5/^m 

8.0 - 8.7/am 

Havana, Cuba 9-29-58 
H 2 0 = 4.25 cm 

299.00 

288.21 

292.52 

290.04 

Kenira, Morocco 9-29-58 
H 2 0 = 2.89 cm 

293.00 

286.79 

290.19 

287.32 

Oakland, California 9-29-58 

H O = 2.60 cm 
2 

289.00 

288.65 

290.00 

287.37 

Maniwaki, Quebec 9-29-58 

HO = 1.39 cm 
2 

274.00 

272.35 

273.56 

271.65 

Ft. Smith, Arkansas 9-29-58 
H 2 0 = 0.63 cm 

267.00 

266.73 

266.98 

265.83 

Chaguaramas Bay, Trinidad 4-1-58 

H O = 3.90 cm 
2 

298.00 

288.17 

292.15 

289.91 

Washington, D. C., 4-1-58 

H O = 0.97 cm 
2 

278.00 

275.91 

277.31 

275.24 

Kindley Berm, Bermuda 8-1-58 
H 2 0 = 5.08 cm 

300.00 

287.58 

291.93 

289.88 

Shreveport, Louisiana, 8-1-58 
H 2 0 = 3.79 cm 

298.00 

291.40 

294.87 

292.16 

Portland, Me. 8-1-58 
H 2 0 = 1.78 cm 

292.00 

288.45 

290.76 

288.15 

Eureka, N.W.T. . 8-1-58 
H 2 0 = 1.28 cm 

278.00 

275.74 

277.20 

275.26 

.[ 





FIGURE CAPTIONS 


Figure la. This illustrates the variation of the difference between window 
channel measurements as a function of cloud amount, if clouds 
radiate as black bodies. 

lb. The effect of decreasing the emissivity at 3.7/xm is shown in 
figure lb. 

Figure 2. Calculated values of T 3 - T n as a function of cloud amount, for a 
cloud emissivity of .75 at 3.7/xm. With a high percentage of 270° 
cloud cover it is not possible to distinguish between T s = 280° and 
T s = 290°. 
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Figure la. This illustrates the verification of the difference 
between window channel measurements as a function of 
cloud amount, if clouds radiate as black bodies. 



Figure lb. The effect of decreasing the emissivity at 
3.7/^m is shown in figure lb. 
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